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O 

, Using a microscopic phase-space model ol the membrane system, the boundary condition at 

^s^j ■ a membrane is derived. According to the condition, the substance flow across the membrane is 

proportional to the difference of the substance concentrations at the opposite membrane surfaces. 
The Green's function of the diffusion equation is found for the derived boundary condition and the 
exact solution of the equation is given. 

PACS number(s): 82.65.Fr, 66.10.Cb 

I. INTRODUCTION 

The membrane systems play an important role in several fields of technology |l]], where the membranes are used 
as filters, and biophysics where the membrane transport is crucial for the cell physiology. The diffusion in the 
membrane system is also interesting by itself as a nontrivial stochastic problem. The process is usually described by 
means of the diffusion equation. Since the equation is of the second order one needs two boundary conditions at the 
membrane surfaces. One condition is provided by the particle number conservation but there is no obvious choice of 
the second condition. Although several proposals can be found in the literature ||,|| a systematic derivation of the 
d ' boundary condition within a microscopic model of the process is missing. 

While the problem of diffusion in a system with a partially permeable wall (the membrane) is not well explored, 
there is a vast literature devoted to the systems with (partially) reflecting or (partially) absorbing walls. Two classes 
of microscopic approaches have been repeatedly used to derive the boundary condition for these cases: the random 
walk on the discrete lattice, see e.g. I^-Q, and the diffusion in the (continuous) phase-space, see e.g. |||||. The first 
approach starts with the death-birth equation while the second one usually utilizes the Fokker-Planck equation. 

In this note we refer to a simple phase-space model of the membrane system and, following P||, we derive, for the 
first time to our best knowledge, the boundary condition for the diffusion equation. The condition appears to be rather 
natural - the substance flow across the membrane is proportional to the difference of the substance concentrations at 
the membrane. In the second part of our paper we find the Green's function of the diffusion equation which satisfy 
the derived boundary condition. The exact solution of the equation with a specific initial condition is also given. 
Finally, our considerations are critically discussed. 



II. DERIVATION OF BOUNDARY CONDITION 



Let consider a system where the substance diffuses across the membrane which is treated as an infinitely thin 
partially permeable wall located at x = 0. Since the concentration gradients are assumed to be only along the 
x— direction the problem is effectively one dimensional. One can imagine the membrane as a reflecting wall with 
homogeneously distributed holes. Then, the substance particles are either reflected by the wall or pass through it 
without change of their momenta. We denote by a the ratio of the total surface of all holes to the membrane surface. 
Then, the probability to pass through the membrane just equals a. The limiting case a = corresponds to the 
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holeless, not permeable wall while a = 1 to the fully permeable wall or to lack of membrane at all. The parameter a 
is, obviously, independent within our model of the particle velocity . 

The distribution function of the substance particles f(p, x, t), where p and x are, respectively, the particle momentum 
and position at the time t, satisfies the relations 

f(-p,0-,t) = (l-a) f(p,0-,t)+af(-p,0+,t) , (1) 

f(p, 0+ t) = (1 - a) f(-p, 0+ t) + a f(p, 0- t) . (2) 

Defining the particle flow and the partial flows as 

j(x,t) =' / dp — f(p,x,t) 

J-oo m 

j+(x, t) d = / dp — f(p, x, t) (3) 

Jo m 

3-(x,t) = ~ [ dp£f(p,x,t) (4) 

J-oo m 

with m being the particle mass, one rewrites the relations (QJ^) as 

j-(0-,t) = {l-a)j+(0-,t)+aj-(0 + ,t) , (5) 
i+(0+, *) = (1 - a) j-(0+,t) + a j+(0- t) . (6) 



One observes that adding Eqs. (|5|,|6|) and using the formula j(x, t) = j+(x, t) — j-(x, t) we get the conservation of the 
substance flow at the membrane i.e. 

j(0-,t)=i(0+,t). (7) 
Now, we decompose the distribution function as 

f( P ,x,t)=f (p) + f 1 (p,x,t), (8) 
where fo(p) is equilibrium distribution function i.e. 

n(x,t) _ \ p 2 



fo(p,x,t) = ' exp 



with T and ks denoting the temperature and Boltzmann constant, respectively. In contrast to the equilibrium 
distribution function, which is even (/o(— p) = fo(p)), the function /i is assumed to be odd (fi(—p, x, t) = —fi(p, x, t)). 
The decomposition (g), which plays a crucial role in our considerations, can be justified in two closely related ways. 
One can assume that f(p,x,t) satisfies the Fokker-Planck equation. Then, the decomposition (||) corresponds to the 
first two terms of the expansion of the distribution function in the large friction limit || . One can also refer here to 
a more general scheme of Chapman-Enskog expansion which is applicable not only for the Fokker-Planck equation. 
However, the function fx is not obligatory odd in this case. Nevertheless, this is the odd part of /i which really 
matters; the even part does not contribute to the particle flow while its contribution to the particle density can be 
neglected because fo ^> \fi\- At the end of our paper we briefly discuss how our results would be modified if the 
decomposition (||) is not limited to the two terms. 

Substituting the distribution function of the form (|^) into Eqs. and applying the Fick law 

j(x,t) = -D d -^l, (9) 

where D is the diffusion constant, one finds 



=0 = 



: - _ / k sT i ^ ( D dn(x,t) 



2irm ' 2 dx 



i-(0 ± 1 t) = A/l£-n(0 ± ,t) + -^Z (ii) 



x=0 ± 
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The relations ©0) allow one to convert Eqs. (0,0) to the form 



,^g(n(r lt )-^ lt )) 

,/g(,(0 +) *)-n(0-,*))= + 



2 — a dn(x, t) 



dx 

2 — a dn(x, t) 



a D dn(x,t) 
x=o- 2 dx 

a dn(x, t) 
x=o+ 2 dx 



x=0+ 



x=0- 



2 dx 

As seen, Eqs. (|12J,|13|) provide again the flow conservation (Q) and the desired boundary condition 

j(0,t) = - K (n(0+ t)-n(0-,*)) 

where 

i(o,t)=i(o+ t)=j(o-,t) , 

and the membrane permeability coefficients /c is defined as 



a / fcsT 



1 — a V 27rm 



III. SOLVING THE DIFFUSION EQUATION 

In this section we are going to solve the diffusion equation 

dn(x, t) d 2 n(x, t) 
dt dx 2 

with the boundary conditions and the initial one given as 

n(x, 0) = no(x) . 

We introduce, as usually, the Green's function G(x,x' ,t,t') which solves the equation 

d , ) d 2 

— G(x, x , t, t) = D — G(x, x , t, t) , 



with the initial condition 



and four boundary conditions: two which are obtained from (uLMJ) i.e 



lim G(x,x',t,t') = S(x-x') 



— G(x,x',t,t') = — G(x,x',t,t') 
dx x=o- dx 



d 



x=0~ 



and 



D— G(x,x',t, t') = -k \G(0 + ,t;x',t') - G(0~ ,t;x' ,t')] , 
dx x=o- 



lim G(x,x',t,t') = 

x — >±oo 



Having the Green's function the solution of the equation (|T^) is given as 

n(x,t) — J dx' no(x') G(x,x',t,Q). 
Further, we always put t' = and denote the Green's function as G(x, x' ,t). 



We find the Green's function applying the standard procedure pOl of the Laplace transformation which gives 



G{x, x', s) = 



dt , 



G(x,x',t) 



Transforming Eq. (h_7|) and the initial condition (h8|) we get the equation 



D— G(x, x', s)-s G(x, x', s) = -5(x - x') . (23) 
dx z 

One trivially finds the general solution of the homogenous equation (^3|) while the inhomogenous equation is easily 
solved by means of the Fourier transformation. In this way we get 



G(x, x , s) — 



1 



2Dq 



- q \x-x'\ + Ae qx + B. 



where q 2 = s/D while the constants A and B are determined by the boundary conditions ( |l9| , |20| , 21). To calculate 
the constants one has to distinguish four cases related to the signs of x and x' . We denote as G ++ the Green's function 

which corresponds to x > and x' > 0; as G-| that one for the case of x > and x 1 < 0, etc. After inverting the 

Laplace transform we get 



G A 



.(x,x',t) = G-+(x,x',t) = — exp(^ 



k /2k(\x - x'\ +2nt) 



D 



erfc 



\x - x'\ + 4Kt 



G++(x. x' , t) — G--(x, x' , t) 



1 



exp 



/ {x-x'f \ 
\ ADt ) 



exp 



2VDt 
(x + x') 2 



k /2k(\x + x'\ + 2nt) 
D eXP ( 



D 



^ erfc^ 



ADt 

\x + x'\ + AkI 



(24) 
(25) 



2VDt 



with erfc(a;) being the complementary error function i.e. 



erfc (it) 



dt e 



As an application of the Green's functions (24 2^) we consider the time evolution of the concentration of the 
substance which is initially homogeneously distributed in the left half-space i.e. the initial condition reads 



n(x,0) — no 0(— x) . 



Then, Eq. ( f22| ) provides 



n + (x,t)—no / dx' Ga (x,x',t) 



2 



CrfC (OT)~ eXP ( 



(2n{x + 2nt) 



D 



erfc 



' x + Ant\ 
> 2-JWt ' 



(26) 



n_ (x, t) = no / dx' G__ (x, x' , t) 



no 
2 



2 — erfcf — ^= ] + exp( 

\2Vd1J V 



2k(-x + 2nt) 
D 



erfc^ 



-X + Ant 
2-JlTt 



(27) 



where n + (x,t) = n{x,t) for x > and n {x, i) = n(x,t) for x < 0. 

Let us briefly discuss the solution (|2(], |27]). One observes that when k = 0, which corresponds to the fully reflecting 
wall, we get from ( p6| , |27j ) an expected result i.e. n + {x,t) = and n_(x,t) = n a . In the opposite limit k — > oo 
corresponding to the lack of any membrane, we find the well known formulas: 
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n+(x, t) = ^ erfc^^— ^== ) 



n-{x,t) = — 



Finally, one observes that for \f~Dt ~^>\x\ and nt^> \x\ the solution (p^, [2^) provides 



n±(x,t) 



"o 



it 



1 nD 
2k V ~T~ 



Thus, n±(x,t) = uq/2 for t — > oo. 



IV. DISCUSSION 



The decomposition (g|) , which, as already mentioned, plays a key role in the boundary condition derivation, includes 
only two terms. Since the solution of the Fokker-Planck equation can be systematically expanded in the powers of 
the inverse friction coefficient, one can easily supplement, following 0, the decomposition (g) by a next order term. 
Assuming that the Fick law (Q) still holds, we again get the current conservation (Q) while the equation analogous to 
(H) reads: 

a D ( d 2 n(x,t) d 2 n(xA) \ / , , . , .\ , , 



where 7 is the friction coefficient from the Fokker-Planck equation. If, following j8|, one treats LHS of Eq. (28) as 
first two terms of the Taylor expansion of the current at a finite x, the boundary condition ( p8[ ) can be manipulated 
to the form 

(1 " fc) ^"T^' i ) +A:j (T' t ) = ~ K ( n(0+ ' i)_n(0_ ' i) ) ' (29) 

where k is an arbitrary number while xq denotes the thickness of the kinetic near membrane layer equal 

a 

x = 



4y / 7r (1 — a) 7 



Unfortunately, we can not see a way to remove the arbitrarness of k which appears when, due to Eq. (Q), we 
express j(0,t) as (1 — k)j(0 + ,t) + kj(0~,t). Such a problem is absent in the formula analogous to ( |29| ) for the case 
of absorbing wall ||. The results (^8[ |2^) are interesting by themselves but they do not seem to be very useful as 
membrane boundary conditions. As well known, the diffusion equation provides a reliable description if the Chapman- 
Enskog expansion converges fast. Then, the third term is not needed in the decomposition (||) and one should use the 
boundary condition (|lj). However, the distribution function can significantly deviate from the equilibrium form in 
the near membrane layer as it happens in the vicinity of the absorbing wall [jll| . Then, the boundary condition ( |29| ) 
seems to be natural but the diffusion equation should be combined with the Fokker-Plnack for a correct description 
of the kinetic layer and the whole approach appears to be very cumbersome. 

The formula of the membrane permeability coefficient ( |l5|) suggests that the permeability grows proportionally to 
\/T, However, one should remember that this dependence has been obtained within a highly simplified mechanical 
model of the membrane which is treated as an infinitely thin reflecting wall with homogeneously distributed holes. 
Then, the substance particles are either reflected by the wall or pass through it without change of their momenta. 
In reality, the membrane structure is much more complicated and the parameter a should be treated as an effective 
probability to pass through the membrane. Consequently, a might be temperature dependent and then k is no longer 
proportional to \/T. 

In the series of papers of one of us j|, the boundary condition different than ( |l4| ) has been advocated. Namely, 
using the symmetry arguments of the Green's function, there has been found instead of ( |l4| ) the relation: 

j(0,t) = (l-<J)io(0,t), (30) 

where S is a dimensionless membrane permeability coefficient and jo is the substance flow in the system with the 
removed membrane. As discussed in [Q], there is a transparent probabilistic interpretation of the relation (]3(]). The 
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solution of the diffusion equation, which satisfies the condition (30) |4|], is qualitatively different than that found here. 
According to the solution, there is a finite flow across the membrane at infinite time while Eqs. ( ^6| , p7|) give the 
vanishing flow in this limit. 

We hope that an experiment will help to chose a right boundary condition and we plan to perform such an analysis 
in collaboration with experimentalists. However, one should remember that there is a whole variety of the mechanisms 
of the substance transport across the membrane jy . If the boundary condition depends on the actual mechanism the 
problem does not have a unique solution. 
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